Surface Outflow in Optically Thick Dust Disks by the Radiation 

Pressure 1 



O 
o 

(N 
u 

< 

(N 



> 
O 

in 

o 
o 

6 



13 



Taku Takeuchi and D. N. C. Lin 

XJ CO /Lick Observatory, University of California, Santa Cruz, CA95064 
taku@ucolick . org, lin@ucolick . org 



ABSTRACT 

We study the outflow of dust particles on the surface layers of optically thick disks. At 
the surface of disks around young stars, small dust particles (size < 10 /jui) experience stellar 
radiation pressure support and orbit more slowly than the surrounding gas. The resulting tail- 
wind imparts energy and angular momentum to the dust particles, moving them outward. This 
outflow occurs in the thin surface layer of the disk that is exposed to starlight, and the outward 
mass flux is carried primarily by particles of size ~ 0.1 /im. Beneath the irradiated surface 
layer, dust particles experience a head-wind, which drives them inward. For the specific case of a 
minimum-mass-solar-nebula, less than a thousandth of the dust mass experiences outward flow. 
If the stellar luminosity is 15 times brighter than the sun, however, or if the gas disk mass is 
as small as ~ 100M®, then the surface outflow can dominate the inward flux in certain radial 
ranges, leading to the formation of rings or gaps in the dust disks. 

Subject headings: accretion, accretion disks — planetary systems: formation — solar system: formation 



1. Introduction 

Circumstellar disks around young stars are 
composed of gas and dust particles, and in gen- 
eral, the dust and gas accretion velocities are not 
equal. The gas suffers turbulent viscosity and ac- 
cretes onto the central star on a viscous time-scale 
of order 10 6 - 10 7 yr (Hartmann et al. 1998). The 
orbital evolution of dust particles, on the other 
hand, is controlled by gas drag (Adachi, Hayashi, 
k Nakazawa 1976; Weidenschilling 1977). 

Several factors affect the radial motion of parti- 
cles. First, gas drag in the radial direction entrains 
particles in the accretion flow. When the gas den- 
sity is high or when the particle size is small, the 
particles follow the gas motion, and accrete onto 
the star on the viscous time-scale. Second, angu- 
lar momentum is transfered between particles and 
the gas. The radial gas pressure gradient shifts 
the angular velocity of the gas slightly from the 
Keplerian value. Moreover, if starlight can ir- 
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radiate particles, then radiation pressure causes 
non-Keplerian deviations in the particle velocity 
as well. Dust particles can thus gain or lose an- 
gular momentum through gas drag and thus move 
outward or inward, depending on whether they ex- 
perience a net head-wind or a net tail-wind. 

In a previous paper (Takeuchi & Lin 2002, here- 
after Paper I) , we studied particle motion in opti- 
cally thick disks, using the assumption that ra- 
diation pressure on the dust is negligible. We 
found that dust particles smaller than 10 /im ac- 
crete with smaller velocity than the gas. As the 
gas accretion proceeds, small dust particles are left 
behind in the disk and the dust-to-gas ratio gradu- 
ally increases. Because circumstellar disks associ- 
ated with young stars are generally optically thick, 
the majority of dust particles are not directly ex- 
posed to the stellar radiation, and the approxi- 
mation used in Paper I is applicable. There are 
always, however, some particles near the disk sur- 
faces that are directly irradiated, and the motion 
of such particles are affected by radiation pressure. 
Klahr & Lin (2001) and Takeuchi & Artymow- 
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icz (2001) studied the motion of particles in opti- 
cally thin gas disks around young Vega-type stars, 
where the dust experiences both gas drag and ra- 
diation pressure. Because the radiation pressure 
reduces the angular velocity of the dust below that 
of the gas, the dust tends to move outward. Thus, 
in a more refined description of an optically thick 
disk, one expects that particles in the irradiated 
surface layer move outward, while beneath the sur- 
face layer particles move inward. If the surface 
outflow dominates the midplanc inflow, then there 
is a net radial expansion that may induce the for- 
mation of rings or gaps. 

In this paper, we study the surface dust outflow 
in optically thick disks due to radiation pressure. 
The mass flux of the surface outflow is calculated, 
which is compared to that of the inflow inside the 
disk. We find that the surface outflow is negligible 
in a minimum-mass-solar-nebula model, but that 
the surface outflow can dominate inflow at some 
disk radii in some cases, including situations in- 
volving high stellar luminosity or low gas mass. 
In such disks, the appearance of rings or gaps is 
expected. 

2. The Distribution of the Gas and the 
Dust 

2.1. The Gas Disk 

We consider circumstellar disks around young 
pre-main-sequence stars. The disks are turbulent 
and are accreting onto the central stars, but they 
are assumed to be passive, that is, the disks are 
primarily heated by stellar radiation only. When 
turbulent fluctuations are averaged out, the mean 
gas motion is nearly in hydrostatic equilibrium 
and the radial accretion flow is much slower than 
the azimuthal Kepler rotation. In cylindrical co- 
ordinates (r, 9, z), the balance of forces in the r- 
and z-directions are 
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where fl g is the gas angular velocity, M is the cen- 
tral star's mass, p g is the gas density, and P g is the 
gas pressure. For simplicity, we assume that the 
gas disk has power-law density and temperature 



profiles in the radial direction, and is isothermal 
in the vertical direction. The temperature distri- 
bution is written as a power-law of index q, 



T(r, z) = T r q Al! 



(3) 



where the subscript "0" denotes quantities at 1 
AU, and a non-dimensional quantity tau is the 
radius in AU. The isothermal sound speed is c = 
c o r Au an d the gas disk scale height is thus defined 
as 



h g (r) = 



(4) 



^K,mid 

where £l K = [GM/(r 2 + z 2 ) 3 / 2 ] 1 / 2 is the Keplerian 
angular velocity and the subscript "mid" indicates 
the midplane value. From the balance between 
the stellar gravity and the gas pressure gradient 
in the z-direction (eq. [2]), we have the density 
distribution 



p g (r, z) = p r p AV cxp (-7^2) 



(5) 



where a power-law in r is assumed. The surface 
density of the gas disk is given by 



S„(r) 



+ OO 



p g dz = V2Trp h Q r A s v , (6) 



where p s = p + (q + 3) /2. The angular velocity Q g 
of the gas is slightly different from the Keplerian 
angular velocity Ok because of the radial pressure 
gradient in the gas. From equation (1), we have 



n 2 (r,z) = (l-r,)ni(r,z) 



(7) 



where f] — —(rfl 2 l p g )~ 1 dP g /dr is the ratio of the 
gas pressure gradient to the stellar gravity in the 
radial direction. As derived in equation (17) in 
Paper I, r\ is written to order (h g /r) 2 as 



v = - 



p + q + 



q + 3 z 2 
2 hi 



(8) 



The disk has a turbulent viscosity v. Turbulent 
disk viscosity has been ascribed to various mecha- 
nisms, including the convective instability (Lin & 
Papaloizou 1980) and the magneto-rotational in- 
stability (Balbus & Hawley 1991). In this paper, 
we simply model the viscous effect of turbulence 
using the so-called a. prescription (Shakura & Sun- 
yaev 1973), 



v = ach g = ac h a r A ^^ 2 



(9) 
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We adopt the following fiducial parameters: M — 
1 Mq, po = 2.83 x 10- 10 g cm- 3 , h a = 3.33 x f CT 2 
AU, p = -2.25, q = -0.5, and a = 10~ 3 . These 
values correspond to a disk having a gas mass of 
2.5 x 1CT 2 M Q inside 100 AU. The power law index 
of the surface density is p s = p+(q + '3)/2 = —1.0. 

— 1/2 

The temperature distribution is T = 278 r AU ' K. 

2.2. The Dust Disk 

2.2.1. Radial Distribution 

The radial surface density distribution of the 
dust is assumed to be proportional to that of the 
gas disk, so that E d)a ii(r) = /duster), where 
the constant /dust is the dust-to-gas mass ratio. 
Throughout the paper, we adopt the value /dust = 
0.01 at all radii. The subscript "all" represents 
contributions from the entire distribution of par- 
ticle sizes. In this paper, we consider cases where 
all particles have the same size and also power- 
law size distributions. For power-law size distribu- 
tions, the size range is bracketed at the large and 
small ends, [s m ; n , s max ]. We adopt a power-law in- 
dex —3.5, appropriate to the distribution expected 
following collisional evolution (Hellyer 1970). The 
surface number density of particles in a size range 
[s, s + ds] at a distance r from the star is given by 

N? : (r,s)ds = N 1 (r)s-Z- 5 ds , (10) 

where s cm is the particle size in centimeters. The 
radial distribution of particles N\ is 

N i(r) = ^-/dustSg(r) (s^ cm - s^i 2 n cm ) (cm)" 4 , 

(11) 

where p p is the physical density of a particle. We 
adopt p p = 1.25 g cm" 3 , / dus t = 0.01, s min = 
0.01 pm, and s max = 1 pm. 

2.2.2. Vertical Distribution 

Dust particles feel the z-component of the stel- 
lar gravity, and thus sediment toward the disk 
midplane. Small particles (s < 1 cm) immediately 
(within one orbit) reach the terminal velocity v z ^, 
where the gravity and the gas drag are in balance. 
For small particles considered here (s < 1 cm), the 
Epstein's drag law is applied, because the mean 
free path of the gas molecules (7 cm at 1 AU) is 
larger than the particle size. The gas drag force is 

F g = -^np g s 2 v t (v d - v g ) , (12) 



where Vt is the mean thermal velocity of the gas 
molecules, and Vd and v g are the velocities of the 
dust and the gas, respectively. The stopping time 
due to the gas drag is t s = md\vd—v g \/\F g \, where 
tnd = |ttp p s 3 is the particle mass. We define the 
non-dimensional stopping time normalized by the 
orbital time as, 

rp , Pp Stl K,mid /-.„n 

I s — ts^K.mid — , U-Jj 

PgTVt 

where «K,mid is the Keplerian velocity at the mid- 
plane. When terminal velocity is achieved, the 
gas drag balances the z-componcnt of the stel- 
lar gravity, —md^l^z. Using the approximation 
12k ~ ^K,mid, the terminal velocity v z .d is 

V Z A = -^K,mid?> • (14) 

The mass flux of particles in the vertical direction 
is 

/ Z ,scd = PdV z ,d , (15) 

where pd is the dust mass density. If the gas disk 
were laminar, the particles would concentrate at 
the midplane in a time that is short compared to 
the evolution time of the gas disk (the sedimenta- 
tion time-scale is about 10 6 yr for 1 pm. particles 
at 10 AU). When the gas disk is turbulent, how- 
ever, the gas stirs and disperses the particles and 
prevents sedimentation. The mass flux of dust 
arising from this turbulent diffusion is assumed 
proportional to the gradient of the dust concen- 
tration pd I p g in analogy with molecular diffusion 
(e.g, Monin & Yaglom, 1971, p. 579; Morrill, 1985, 
§3.5.1; Cuzzi, Dobrovolskis, & Champney 1993). 
In the vertical direction, the mass flux / Zj dif is 

, p g v d f p d \ 

fzdit = -^d-z{7g) ' (16) 

where the Schmidt number, Sc, represents the cou- 
pling strength between the particles and the gas. 
For small particles that are well coupled to the 
gas, Sc approaches unity. For large particles that 
are decoupled from the gas, Sc becomes infinite. 
In our standard model, we use Sc = 1. 

In equilibrium, the sedimentation mass flux bal- 
ances the diffusive mass flux: 

/ z ,scd + /,,dif = . (17) 

The equilibrium dust density distribution is given 
by equation (31) in Paper I. For particles in a size 
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range [s, s + ds], the density distribution is 

•v2 



Pd(r,z,s)ds 



Pd(r, 0, s) exp 



exp 



(18) 



where X^mid is the midplane value of the non- 
dimensional stopping time. The total dust density 
is s 

Pd. A \\{r,z) = I p d {r 1 z,s)ds . (19) 

As discussed in Paper I, the sedimentation of 
particles is more effective for larger particles and 
at larger distances from the star. As increasing the 
distance, the dust disk becomes relatively thinner 
compared to the gas disk. 

2.3. Optical Depth from the Star 

The disks considered in this paper are optically 
thick at visible wavelengths even in the vertical di- 
rection (t voi . <~ 100 at 100 AU). Stellar light cannot 
penetrate the dust disk. Only particles residing in 
the thin surface layer are exposed to stellar radia- 
tion pressure. 

The optical depth from the star to the position 
(r, z) in the dust disk is 



r(r,z)= f J^ PdMl (r',z')Jl+ Z ^dr' , (20) 

Jr in /dust V H 

where z' = (z/r)r', r in is inner edge of the dust 
disk, and n is the dust opacity in the visible wave- 
length for an unit mass of the gas. The opacity 
for a unit mass of the dust is ^/./dust- In our cal- 
culations, we adopt r; n = 0.1 AU. This choice of 
inner boundary does not affect the integral (20), 
because the dust disk is flared and the main contri- 
bution to the optical depth comes from the dust at 
relatively large distances from the star. We adopt 
a dust mass opacity n = 100 cm 2 g _1 for visible 
wavelengths. For particles smaller than 100 /zm, 
the opacity is independent of (the maximum) par- 
ticle size (see Figs. 5 and 11 in Miyake & Naka- 
gawa 1993). Thus, we need not consider variations 
in k with changes in the particle size. 

Figure 1 shows the disk surfaces at which the 
optical depth to the starlight is unity. Above these 



illuminated surfaces, the dust particles are directly 
exposed to the stellar radiation. Figure la shows 
the illuminated surfaces of the standard disk mod- 
els which are composed of single-sized particles 
(0.1 to 10 fim). As the particle size increases, the 
level of the illuminated surfaces becomes lower. 
This size dependence is because larger particles 
sediment more readily, and thus concentrate near 
the midplane. The illuminated surfaces have a 
slightly flatter geometry than the gas disk, be- 
cause the dust particles at larger distances sedi- 
ment more effectively. Figure 16 shows the illumi- 
nated surfaces of the models containing power-law 
size distributions. The illuminated surfaces lie at 
roughly three disk scale heights, which are similar 
to those of the single-sized 0.1 /im particles. The 
location of the illuminated surface does not change 
significantly with the maximum size of the parti- 
cles, because high altitude above the midplane is 
mostly populated by the smallest particles. 

3. Outflow in the Surface Layer 

Dust particles above the illuminated surface 
suffer radiation pressure from the central star. 
The radiation pressure balances a part of the stel- 
lar gravity, so that the small particles orbit slower 
than the gas. Because the gas drag transfers an- 
gular momentum from the gas to the particles, the 
particles move outward, leading to an outflow in 
the surface layers of dust disks. In this section, we 
estimate the mass fluxes carried by these outflows. 

3.1. Outflow Mass Flux 

The radiation pressure is expressed through the 
ratio (3 to the stellar gravity. The radial compo- 
nent is 

Frad - (3m d rfl 2 K . (21) 

The value of (}, which is calculated with the Mie's 
scattering theory, depends on the particle compo- 
sition, porosity, shape, and size. In this paper, 
we adopt j3 for "young cometary particles," which 
is calculated by Wilck & Mann (1996) for typical 
interplanetary particles in the solar system and 
shown in Figure 2. In general, stellar light causes 
Poynting- Robertson drag on the particles in addi- 
tion to the gas drag induced by radiation pressure. 
However, in the gas disks considered here, the gas 
density is sufficiently large so that gas drag domi- 
nates Poynting- Robertson drag (see the discussion 
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in Takeuchi & Artymowicz 2001, §4.2). Poynting- 
Robertson drag can therefore be neglected. 
The equations of motion of a particle are 



dv r d v 9.d mr> 2 ^K.mid 



(Vr,d - Vr,g) > 
(22) 



«K,mid 



{V6,d - V0,g) , 



(23) 



where v r and ve are the r- and 9- components 
of the velocity, and the subscripts "g" and "d" 
distinguish gas and dust. As in Paper I, we as- 
sume v e , g w ~ WK,mid and d(rve,d)/dt w 
w r ,dd(™K,mid)/^ = w r ^u K ,mid/2. From equation 
(23) we then have 

«M _ = --jTsV r ,d ■ (24) 

In equation (22), dv r ^ d /dt = 0(v 2 d /r) is neglected 
if («r\d/wK,mid) 2 <C 1. We also neglect second or- 
der terms in r\ and (v 8 ^ d — we, 9 )/«K,mid- To this 
level of approximation, v 2 . d v 2 . g + 2vg g (vg d — 

v 0,g) ~ w k(1 ~V) + 2v K ,mid{ve,d - ve, g ), and equa- 
tion (22) becomes 



(f3-r]) v K+ 2v K,mid(ve,d-ve,g)- VK ^ ld ( Vr 



-V r> g) = . 

(25) 



Using equation (24) and vk ~ «K,mid, the particle 
radial velocity is 

v rA r, Z ,s)= T ^[^> K ^ . (26) 

The value of r\ is of order {h g /r) 2 and it is smaller 
than 0.05. Since particles smaller than 10 pm have 
/3 » 77 (see Fig. 2), r\ in equation (26) can be ne- 
glected. Particles larger than ~ 100 pm may have 
(3 smaller than 77. For such large particles, the ra- 
diation pressure is unimportant, even if they are 
above the illuminated surface. In this section, we 
consider only small particles with radii less than 
10 /ith. Such particles have non-dimensional stop- 
ping times T s « 1 (see Fig. 2 in Paper I). The 
radial velocity then reduces to 

Vr4 = V r ^ g + (3T s V K ,mid ■ (27) 

The mass flux of the outflow in the surface layer 
is 

F aul (r) = 2 2-Krp d (r 1 z 1 s)v r ^ d (r,z,s)dzds , 

J S m i n J -2 sur 

(28) 



where z sm is the height of the illuminated sur- 
face. The factor 2 arises from considering both 
the top and bottom surfaces. In the integration 
in z above the illuminated surface, v Tyd varies less 
rapidly than p d , and we may use its value at the 
illuminated surface. The mass flux is then ex- 
pressed as 



F SUI (r)=2wr v r , d (r, z sur , s)S 

J Smin 



d, sm (r,s)ds , 
(29) 

where T, dtSur ds is the surface density of dust parti- 
cles above the illuminated surface in the size range 
[s, s + ds], 

f-OO 

£ d:Sur (r,s) = 2 p d (r,z,s)dz . (30) 

The surface outflow can clear out the dust within 
a distance r from the star on a time-scale t ou t ~ 
M Aust (r)/F sur (r), where M dust (r) is the dust mass 
within r. In our models, the disk mass is dom- 
inated by the outermost part of the disk and 
M dust (r) w 7rr 2 £^ a n(r). The dust clearing (or 
the disk evolution) time-scale is thus 



nr 2 Y, dM i(r) 



(31) 



While at the surface layer particles flow out- 
ward, the particles beneath the surface layer are 
shielded from stellar irradiation. At such loca- 
tions, the gas drag causes a net inward flow of 
particles as shown in Paper I. From equation (26) 
with p = 0, the radial velocity of particles beneath 
the surface layer is 



v r ,d{r,z,s) 
The mass flux is 



T s 1 V r ,g - ??^K,mid 



(32) 



2irrp d (r, z, s)v r , d (r, z, s)dzds . 

(33) 

The time required for the inflow to transport a 
mass equal to the dust mass within r is 



*in(r) 



7rr 2 S rf . al i(r) 

\Fn{r)\ 



(34) 



We define this as the disk evolution time through 
the inflow. This definition of t- m docs not mean the 
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time for carrying the total dust mass outside r to 
the inner part, but it is convenient for comparing 
with tout- Comparison of t ou t and i; n simply means 
comparison of the outflow and inflow fluxes. The 
disk clearing through the surface outflow occurs 
effectively only if the surface flux is stronger than 
the inflow flux, i.e., t ou t < Un- The comparison of 
these two time-scales is discussed in §3.3 and §3.4 
below. 

3.2. Evolution Time of the Disks through 
the Outflow 

We first discuss the disk clearing time, t out , 
through the surface outflow. In this subsection, 
we adopt the assumption that inflow flux near the 
midplanc is negligible and that only the surface 
flow affects the dust structure within the disk. In 
§3.3 and §3.4, we relax this assumption and con- 
sider more realistic scenarios including the effects 
of inflow near the midplane. 

The particles in the surface layer have out- 
flow velocity, v r ^(r, z sur , s). These particles are 
blown from their orbits to larger radii on a time- 
scale tb — r/v r ,d{r, z sm , s). For disks composed of 
single-sized particles, tb is plotted in Figure 3. The 
blow-off time for 0.1 to 10 /zm particles is about 
10 2 yr at a few AU from the star, and is about 10 5 
yr at 100 AU. The blow-off time is much shorter 
than the inward migration time inside the disk, 
which is discussed in §3.3 below, because (5 ^> ?] 
and because the stopping time T s at the illumi- 
nated surface is much longer than at the midplane. 

In our standard model, the dust disk is highly 
optically thick. Even in the vertical direction at 
100AU, the optical depth is T vcr ^ 100, so the layer 
above the illuminated surface is very thin. Figure 
4 shows the ratio of the surface density of the layer 
to that of the whole disk (£<j,sur,aii/£<i,aii)- As the 
distance from the star increases, the mass ratio 
increases and the surface layer becomes thicker. 
However, even at 100 AU, the mass of the surface 
layer is less than 10~ 3 times the dust disk mass. 
In the case of single-sized particles (Fig. 4a) the 
mass of the surface layer is smaller for disks with 
larger particle sizes. As the particles grow larger, 
they sediment more to the midplane, and the sur- 
face layer becomes thinner. On the other hand, for 
the mixed-sized particles (Fig. 46) , the mass of the 
surface layer does not vary significantly with the 
maximum size. Because smaller particles exist at 



higher altitudes above the midplane, the surface 
layer is composed mainly of the smallest particles, 
and the quantity of smallest particles in the sur- 
face layer does not vary with the maximum size. 

Although Figure 3 implies that the outflow ve- 
locity at the surface layer is high, the outflow mass 
flux is not large because of the small mass con- 
tained in the surface layer. The disk evolution 
time ^ut through the surface outflow is shown in 
Figure 5a. The evolution time of the disk com- 
posed of 0.1 /im particles is about 10 7 yr at 1 AU 
and 10 8 yr at 100 AU. Thus, in 10 7 yr, the dust 
outflow mass which has passed 1 AU becomes com- 
parable to the dust disk mass inside 1 AU. If the 
radial motion of the dust disk occurs only through 
surface outflow, the particles in the innermost part 
of the disk will escape from the region inside 1 AU 
in 10 7 yr, and a hole will appear. The hole ex- 
pands, and in 10 s yr, the radius of the hole reaches 
100 AU. There is, however, an accretion flow un- 
der the surface layer. Inhibition of disk clearing 
by the accretion flow is discussed in §3.3 below. 
As the particle size increases, the evolution time 
increases as well. This is because particles with 
larger sizes concentrate more at the midplane and 
the surface layer becomes thinner. The blow-off 
time tb for a larger particle at its surface layer is 
also larger because of the higher local density of 
the gas at the lower surface layer (Fig. 3). If the 
disk is composed of 10 fim particles, it takes about 
10 9 yr for a disk of radius 100 AU to evolve. 

Figure 56 shows the evolution times for disks 
composed of mixed-sized particles. The disk con- 
taining [0.01 ^im- 1 /an] particles evolves in about 
3 x 10 7 yr at 1 AU and in 3 x 10 8 yr at 100 AU. 
The disk evolution occurs mainly through the out- 
flow of 0.1 /xm particles. While the surface layer 
is composed primarily of the smallest (0.01 /im) 
particles, these particles do not receive radiation 
pressure as effectively as particles with 0.1 ^m (see 
Fig. 2). Thus, the surface outflow is mainly car- 
ried by a minor fraction of particles in the 0.1 /im 
range, and the disk evolution time of the mixed- 
sized particles is much longer than that of the disk 
with the single-sized particles of 0.1 ^tm. The evo- 
lution time increases with the maximum particle 
size. The dependence of the maximum size on 
the evolution time, however, is not as significant 
as in the case of single-sized particles. For disks 
composed of mixed-sized particles, particles of size 
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0.1 /zm form the main contribution to the surface 
outflow. The number density of 0.1 /zm particles 
decreases as the maximum particle size increases. 
However, because the thickness of the surface layer 
is determined by the optical depth from the star 
and the opacity comes mainly from 0.1 /jm par- 
ticles, the total mass of 0.1 /im particles in the 
surface layer does not change significantly (Fig. 
46). Hence, the outflow flux maintains a relatively 
constant value. 

3.3. Comparison with the Inward Migra- 
tion Time 

The surface outflow leads to disk clearing only 
if it dominates the under-surface inflow. The disk 
evolution time-scale, t[ n , through the inflow is 
shown as the dashed lines in Figure 5. As dis- 
cussed in Paper I, for sufficiently small particles 
that are well coupled to the gas, the inflow ve- 
locity is similar to the accretion velocity of the 
gas. In the standard model, particles smaller than 
1 ^im have velocities nearly indistinguishable from 
the gas throughout the disk. Particles of 10 fim 
decouple from the gas motion at large distances 
(r > 50 AU) and have a shorter evolution time 
than the gas and smaller particles. The inflow- 
ing mass flux is much larger (10 times or more) 
than the mass flux of the surface outflow even for 
0.1 /itm particles. Therefore, the net mass flux of 
particles is dominated by the inflow beneath the 
surface layer. In the disks composed of mixed- 
sized particles (Fig. 56), the inflow also dominates 
the net flow. The surface outflow has a very minor 
effect on the disk evolution in our standard model, 
regardless of the particle size. 

3.4. Various Models 

In our standard model, the surface outflow is 
negligible for disk evolution compared to the in- 
flow inside the disk. This flow outcome is due to 
the surface layer being very thin such that its mass 
is a small fraction of the entire dust disk mass. 
The surface outflow may be important under some 
conditions where the surface layer is thicker or the 
outflow is much faster than the standard model. 
In this subsection, we consider several models to 
see under which conditions the surface outflow is 
important. 



3.4-1- High Luminosity 

Massive young stars have a larger luminosity- 
mass ratio. For example, 2.5M Q stars have about 
20L Q /M Q for ~ 10 7 yr (DAntona & Mazzitelli 
1994). Under such a high luminosity-mass ratio, 
the dust particles are exposed to a stronger radi- 
ation pressure (/3 is proportional to L/M). Fig- 
ure 6a shows the evolution time of the disk in the 
case where the stellar luminosity-mass ratio is 15 
times larger than the standard model. The out- 
flow velocity is 15 times faster than in the standard 
model, causing the outflow mass flux to dominate 
the inflow at large distances (r > 200 AU). In 
this model, dust particles which are initially out- 
side 200 AU are expelled to larger distances, while 
particles inside 200 AU are accreted onto the star. 
It is expected that a gap around 200 AU forms in 
10 7 yr. 

3-4-2. Small Disk Mass 

If the total dust mass is smaller, the disk opti- 
cal depth is smaller so that the layer above the 
illuminated surface is thicker. Figure 66 shows 
the evolution time of the disks whose masses are 
0.1 and 0.01 times smaller than the standard disk. 
The mass flux of the surface outflow is almost in- 
dependent of the disk mass because the mass of 
the surface layer, which is determined by the op- 
tical thickness, is approximately constant. There- 
fore, the relative importance of the surface out- 
flow, which is proportional to the mass ratio of the 
surface layer to the whole disk, increases with de- 
creasing the total dust mass. The evolution time 
through the inflow (see the dashed lines) at the 
inner part of the disk is similar to that of the 
standard model, while at the outer part it de- 
creases with decreasing disk mass. As the disk 
mass decreases, the gas density at the outer region 
of the disk becomes insufficient to keep dust par- 
ticles moving together with the gas. The inward 
drift velocity of particles relative to the gas be- 
comes faster for disks with smaller gas mass. For 
such drifting particles, the evolution time is pro- 
portional to r 2 +P~«/ 2 (see Paper I, §3.3.3), which 
is constant with r in our models (p — —2.25 and 
q = —0.5). For the disk with one tenth the mass 
as the standard model, the evolution time through 
the surface outflow is still longer than that of the 
inflow. If the disk is 10~ 2 times smaller than the 
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standard model (the gas mass inside 100 AU is 
2.5 x 1O~ 4 M0 = 83M e ), the surface outflow flux 
is comparable to the inflow flux at around 20 AU. 
Thus, there is no net flow at radii in the neighbor- 
hood of 20 AU. The dust particles initially outside 
20 AU move inward until their migration is termi- 
nated at 20 AU. At the same time, the particles 
initially inside 20 AU continue to move inward and 
accrete onto the star. Therefore, the particle mi- 
gration induces accumulation of dust at 20 AU and 
forms a ring there. Gas disks around young stars 
are expected to dissipate in 10 7 yr (Hartmann et 
al. 1998). When the gas mass is reduced to about 
100M e , a ring appears at a radius of a few tens of 
AU. 

3.4-3- High Diffusivity of Dust 

In the standard model, the diffusivity of dust 
particles subject to turbulent motion of the gas is 
assumed to be the same as the diffusivity of the 
gas, i.e., the Schmidt number, Sc, is unity. Ex- 
periments show that Sc can be as small as 0.1 
and dust particles are much more diffusive than 
the gas, if the Stokes number, St, is of order 10~ 2 
(e.g., Fig. 1 in Cuzzi ct al. 1993). The Stokes 
number is considered to be of the same order as 
the non-dimensional stopping time T s . At radii of 
several hundred AU, T s ~ 10~ 2 for 0.1 /im parti- 
cles at the illuminated surface (z ~ 3/i g ). Thus, 
it is expected that such particles will spread to 
higher altitude, where they are easily blown off. 
Figure 6c shows the evolution time for such par- 
ticles with Sc = 0.1. The dust particles are more 
diffusive and are stirred up to higher altitude by 
turbulence than in the case of Sc = 1. The illu- 
minated surface is located higher, where the gas 
density is lower and the outflow velocity is faster. 
The outflow mass flux, however, is still smaller 
than the inflow mass flux within r < 500 AU. In- 
flow under the surface layer dominates even for the 
particles with high diffusivity. 

3.4-4- Particle Growth 

In the above discussion, we mainly used mod- 
els which have a power-law size distribution with a 
maximum size of 1 /um. The dust particles in cir- 
cumstcllar disks may have evolved to sizes larger 
than 1 ^m (Beckwith & Sargent 1991; Mannings 
& Emerson 1994; D'Alessio, Calvet, & Hartmann 
2001). Figure 5b shows the evolution time of the 



disks as the maximum size of the particles is var- 
ied. The evolution time through the surface out- 
flow does not change significantly as discussed in 
§3.1. The evolution time through inflow inside the 
disk, however, rapidly decreases with the growth 
of particles. At the outer part of the disk where the 
gas density is insufficient to keep the largest par- 
ticles moving together with the gas, the inflow ve- 
locity of the particles is proportional to their size. 
For the size distribution with a power-law index 
—3.5, the dust mass is dominated by the largest 
particles, and the inflow mass flux is proportional 
to the size of the largest particle. The evolution 
time through the inflow is as small as 10 5 yr when 
the largest particle is 1 mm, which is much shorter 
than that through the surface outflow. As the par- 
ticles grow larger, the surface outflow becomes less 
important for the overall disk evolution. 

4. Discussion and Summary 

4.1. Progress in the Disk Evolution 

We compared the surface outflow flux to the 
under-surface inflow flux and discussed conditions 
for disk clearing. In those comparisons, we as- 
sumed the dust density profiles described in §2, 
and neglected time variations of the dust density. 
In real disks, however, the dust density varies as 
the surface outflow proceeds. We discuss how den- 
sity changes affect the subsequent evolution. 

As dust particles above the illuminated surface 
are blown off by the stellar radiation, new particles 
continually diffuse into and replenish this region 
from below. The diffusion time over which the 
surface layer is re-populate is 

w Sc(3^ m 9 n _, ^ 1Q 3 r 3 /2 yr 
v a 

where we use the fact that the distance between 
the surface layer and the midplane is about three 
disk scale heights for a = 10 -3 and Sc = 1. The 
diffusion time (10 6 yr at 100 AU) is much shorter 
than the disk evolution time t out . Thus during 
the disk evolution, the dust in the surface layer 
is continuously replenished, and the surface out- 
flow continues until the dust in the midplane is 
depleted. 

In disks composed of mixed-sized particles, 
small (s < 0.1 fim) particles constitute the ma- 
jorities of the surface outflow. The small particles 
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are rapidly removed from the disk, and are con- 
tinuously replenished through the collisional de- 
struction of larger particles. The collisional time 
scale, t co \ = (irs 2 ridVr C i) , where is the particle 
number density and v rc i is the relative velocity of 
particles, is much shorter than the disk evolution 
time. For example, if we assume that the major- 
ity of the dust is composed of 1 cm particles, the 
relative velocity is about 10 cm s _1 at 1 AU (from 
Fig. 3 in Weidenschilling & Cuzzi 1993), and the 
collisional time is t co \ ~ 10 3 yr. If the dust is 
composed of smaller particles, the collisional time 
is smaller. The outcome of the collisions, how- 
ever, is not well understood, and we do not know 
how many small particles are produced through 
the collisional destruction. In this paper, we sim- 
ply assume that the particle size distribution is 
in an equilibrium state. The minimum size, the 
maximum size, and the power-law index of the 
distribution do not vary with time. In a realis- 
tic situation, however, at least the maximum size 
would grow with time, leading to inhibition of the 
disk clearing, as discussed in §3.4.4. Calculations 
that take account of dust growth are needed to 
solve this issue. 

Because the disk evolution time t out is shorter 
at smaller radii, the dust outflow proceeds from 
the inner disk. The particles blown from the in- 
ner disk accumulate in the outer disk. This pro- 
cess increases the dust-to-gas ratio at large radii. 
We can ask whether this increase in the dust- 
to-gas ratio affects the subsequent disk evolution. 
Suppose that the dust particles inside a radius r 
are removed and an inner hole (or a gap) is pro- 
duced. The removed particles pile up in an an- 
nulus just outside the inner hole. The width of 
the annulus is larger than the disk scale height. 
On the diffusion time t lp , the accumulated par- 
ticles diffuse over a few disk scale heights both 
in the vertical and in the radial directions, and 
during a disk evolution time t ou t > t Tp , the parti- 
cles diffuse even further. Thus, the surface den- 
sity of the annulus increases at most by a fac- 
tor <~ 7rr 2 S dia n/ '(27rrS dia ndr) <r/h g ~ 10, where 
dr ^ h g ~ O.lr is the width of the annulus. Be- 
cause the gas density does not change, the dust- 
to-gas ratio of the annulus becomes < 10 times 
larger. The opacity n, which is defined for a unit 
mass of the gas and is proportional to the dust-to- 
gas ratio, also increases by a factor of order < 10 



times. (If the particles diffuse as far as the inner 
hole radius r in a disk evolution time t out , then the 
hole is refilled and dust removal is inhibited. The 
time-scale for the refill is comparable to the gas ac- 
cretion time-scale r 2 /v, which is shown as the dot- 
ted lines in Figure 5. In the standard model, the 
refill effectively inhibits the dust removal. Only in 
the high stellar luminosity or the low mass disk 
cases, the surface outflow can dominate the refill.) 

Figure 7 shows how the increase in the dust-to- 
gas ratio affects the outflow mass flux. In Figure 7, 
the outflow flux in the cases for /dust = 0.1 (i.e., 
the dust-to-gas ratio is 10 times larger than the 
standard model) and for /dust = 10~ 3 (10 times 
smaller dust-to-gas ratio) are shown with the stan- 
dard model. It is seen that the outflow flux de- 
pends only weakly on the dust-to-gas ratio. The 
change in the dust-to-gas ratio has a slight effect 
on the location of the surface layer, but it hardly 
changes the dust mass of the surface layer. The 
thickness of the surface layer is determined so as 
to make the optical depth be unity, and its mass 
content does not depend on the dust-to-gas ra- 
tio. The variation in the outflow flux arises mainly 
from change in the location of the surface layer and 
from change in the angle of the incident starlight 
to the surface layer, and these effects are small. 
Therefore, the accumulation of blown off particles 
increases the dust-to-gas ratio at the disk edge just 
outside the inner hole, but docs not significantly 
decrease the outflow flux there. Even after the in- 
crease in the dust-to-gas ratio in the outer disk, 
the surface outflow proceeds and the disk evolu- 
tion continues with the time-scale plotted in Fig- 
ures 5 and 6. 

4.2. Gap Formation in Herbig Ae/Be 
Disks 

The phenomenon of surface stripping can ef- 
fectively evolve the disk structure around young 
luminous stars whose luminosity-mass ratios are 
larger than ~ 15 times the solar value, as shown in 
§3.4.1. In the outer part (> 200 AU) of such disks, 
the surface outflow dominates the inflow through 
the midplane, causing a net outward flow. On 
the other hand, within 200 AU, the dust flows in- 
ward onto the central star. Around 200 AU, the 
dust density decreases. The outflow or inflow flux 
around 200 AU is strong enough that it can carry 
the total amount of the dust inside 200 AU in 10 7 
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yr. This flux maintains a relatively constant value 
throughout the evolution even after the dust den- 
sity and the dust-to-gas ratio have decreased (see 
Fig. 7) . Consequently, the dust density declines at 
an almost constant rate. The time-scale for dust 
diffusion to transport dust from the midplane to 
replenish the surface layer is about 3 x 10 6 yr at 
200 AU. We expect the dust density at 200 AU 
will decline considerably in 10 7 yr. 

The formation of gaps can be examined by radio 
observations. At 200 AU, the gas surface density 
in the standard model is S g = 1.8 g cm~ 2 , which 
is optically thin in the vertical direction in radio 
continuum region around 1 mm (ri mm ~ 10~ 2 ). 
Therefore, the density decline will be observed as 
a decrease in the optical depth. 

The disk is highly optically thick in the visible 
or near infrared light (r vis - 200 at 200 AU). Thus, 
the density decrease is difficult to examine directly 
using visible or near-infrared observations. If the 
gap region is shadowed by the inner disk, or if the 
edge of the outer disk is illuminated by the cen- 
tral star, the variation in the surface irradiation 
will probably be observed. There are a few Herbig 
Ae/Be stars that show gaps in the scattered light 
images of their disks. The disk around HD 163296 
has a dark lane at 325 AU (Grady et al. 2000). 
HD 100546 also has a disk with a gap at 250 AU 
(Grady et al. 2001). In order to make direct com- 
parisons with these observations, we need further 
modeling that takes account of the time evolution 
of the dust density and that incorporates simu- 
lated scattered light images. 

In this paper, we discussed gap formation 
through the surface stripping caused by radiation 
pressure. It is useful to discuss how to distinguish 
this mechanism from the other gap formation pro- 
cesses such as the planetary perturbations (Lin & 
Papaloizou 1986; Takeuchi, Miyama & Lin 1996). 
Obviously, the surface stripping itself cannot gen- 
erate non-axisymmetric features, whereas an em- 
bedded planet can. In gas-free disks, a planet 
forms characteristic disk structures that arises 
from orbital resonances (Roques et al. 1994; Der- 
mott et al. 1994; Liou & Zook 1999; Wyatt et 
al. 1999; Ozernoy et al. 2000; Wilner et al. 2002; 
Quillen & Thorndike 2002), and in gas disks the 
planet excites spiral density waves (Miki 1982; 
Sekiya, Miyama, & Hayashi 1987; Bryden et al. 
1999; Kley 1999; Lubow, Seibert, & Artymowicz 



1999; Miyoshi et al. 1999; Tanigawa & Watanabe 
2002). The detection of such non-axisymmetric 
features would rule out surface stripping as a gap 
formation mechanism. The other notable feature 
of surface stripping is that it does not affect the 
gas density structure. The gas disks still have 
monotonic density structure even after the dust 
disks have formed gaps. Therefore, detection of 
gaps in gas disks also rules out surface stripping. 
In order to explain the cause of the gaps around 
HD 163296 and HD 100546, it is important to 
determine with future observations whether their 
gas disks have gaps or not. 

4.3. Summary 

Dust particles that are directly irradiated by 
the star move outward under the combined action 
of gas drag and radiation pressure. In this paper, 
the surface outflow of particles in optically thick 
disks is studied. Our results are as follows. 

1. Particles in the surface layer, which are di- 
rectly exposed to the starlight, move outward. In 
a disk resembling the minimum-mass-solar-nebula, 
the height of the surface layer is about three scale 
heights of the gas disk and the mass is 10~ 6 times 
the total dust mass at 1 AU and 10~ 4 times the 
total dust mass at 100 AU. For this case, the mass 
flux of the surface outflow is much smaller than 
the inflow mass flux inside the disk, and the sur- 
face outflow does not contribute significantly to 
the dust disk evolution. 

2. The surface outflow dominates the inflow 
in the following cases: (a) If the luminosity-mass 
ratio of the central star is 15 times larger than the 
solar value, the velocity of the surface outflow is 
large enough at the outer part of the disk (r > 200 
AU). (b) If the mass of the gas disk is as small 
as 10 -2 times the minimum-mass disk (~ 100M® 
inside 100 AU), the mass ratio of the surface layer 
to the whole dust disk is large enough at a few 
tens AU. In such cases, the surface outflow induces 
formation of gaps or rings. 

3. The surface outflow is mainly composed of 
particles with size of order ~ 0.1 ^im. As the par- 
ticles grow, the largest particles around the mid- 
plane move inward more rapidly, while the out- 
flow flux at the surface remains at a near-constant 
value. This results in a diminishing importance of 
the surface outflow for the overall disk evolution. 
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Fig. 2. — Radiation pressure to gravity ratio (3 of 
the "young comctary particles." The value is from 
Wilck & Mann (1996). 
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Fig. 1. — Locations of the illuminated surfaces, 
(a) The dust particles are assumed to be single- 
sized. The sizes are 0.1, 1, and 10 ^m from the 
upper line. The dotted line shows the scale height 
of the gas disk, (b) The size distribution is a 
power-law with an index —3.5. The minimum size 
is 0.01 /itm. The maximum sizes are 1, 10, and 
100 /im from the upper line. 
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Fig. 3. — Blow-off time of the outflowing particles 
in the surface layer. The particles are single-sized 
of 10, 1, and 0.1 /im from the upper line. 
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(b) power low size distribution 
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Fig. 4. — Mass ratio of the surface layer above 
the illuminated surface to the whole dust disk 
(£d,sur,aii/£d,aii)- (a) The dust particles are single- 
sized. The sizes are 0.1, 1, and 10 /xm from the 
upper line, (b) The size distribution is a power-law 
with an index —3.5. The minimum size is 0.01 /im. 
The maximum sizes are 1, 10, and 100 /im from 
the upper line. 
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Fig. 5. — Disk evolution time through the outflow 
and the inflow. The solid lines show the evolution 
time through the surface outflow, and the dashed 
lines show the evolution time through the inflow 
inside the disk. The dotted line shows the accre- 
tion time-scale of the gas. (a) The particles are 
single-sized. The sizes are 10, 1, and 0.1 ^im from 
the upper lines (at 10 AU for the dashed lines) . (b) 
The size distribution is a power-law with an index 
—3.5 and a minimum size 0.01 /im. The maximum 
sizes are 1 mm, 100, 10, and 1 /im from the upper 
solid line and from the lower dashed line. 
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Fig. 6. — Disk evolution time for various mod- 
els. The solid and dashed lines correspond to the 
evolution time through the surface outflow and 
through the inflow inside the disk, respectively. 
Particles have a power-law size distribution with 
Smax = 1 M m - Parameters are the same as the 
standard model except for: (a) The luminosity- 
mass ratio of the central star is L/M = 15L©/M©. 
(6) The disk is less massive. The gas masses inside 
100 AU are 2.5 x 10~ 3 M Q and 2.5 x 10- 4 M Q = 
83M ffi from the upper line, (c) The particles is 
more diffusive (Sc = 0.1). 
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Fig. 6 — continued 
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Fig. 7. — Surface dust flux for various dust-to-gas 
ratios, /dust = 0.1, 0.01, and 10 -3 from the upper 
line. The flux is normalized by the value of the 
standard model (/dust = 0.01) at 10 AU. 
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